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Abstract 

The four-year oscillations of the number of spawning sock- 
eye salmon ( Oncorhynchus nerka) that return to their na- 
tive stream within the Fraser River basin in Canada are 
a striking example of population oscillations. The period 
of the oscillation corresponds to the dominant generation 
time of these fish. Various - not fully convincing - ex- 
planations for these oscillations have been proposed, in- 
cluding stochastic influences, depensatory flshing, or ge- 
netic effects. Here, we show that the oscillations can be 
explained as a stable dynamical attractor of the popula- 
tion dynamics, resulting from a strong resonance near a 
Neimark Sacker bifurcation. This explains not only the 
long-term persistence of these oscillations, but also repro- 
duces correctly the empirical sequence of salmon abun- 
dance within one period of the oscillations. Furthermore, 
it explains the observation that these oscillations occur 
only in sockeye stocks originating from large oligotrophic 
lakes, and that they are usually not observed in salmon 
species that have a longer generation time. 

Keywords: cyclic dominance - population oscillations 
- non-linear dynamics - Neimark Sacker bifurcation - 
strong resonance 

1 Introduction 

In many ecological systems distinct population oscillations 
are known, with their specific dynamical characteristics 
often captured by simple generic models. Among these 
are the spatial synchronisation of the lynx-hare oscilla- 
tions in Canada l^f^lE the chaotic oscillations of boreal 
rodents in Fennoscandia'^, or the cyclic outbreak dynamics 
of the spruce budwormf^EI, xhe four-year oscillations of 
sockeye salmon ( Oncorhynchus nerka) in the Fraser River 
basin in Canada are another well-documented example 
of such large-scale population oscillations IH HI Every 
fourth year the abundance of these fish is at very high 
levels, reaching several million fish in some spawning pop- 
ulations, but drops to numbers between several hundred 
and a few ten thousand individuals in the following years 



(hence the term cyclic dominance). Different stocks can 
have their population maximum in different years (Fig. 1 
and Fig. Al in Appendix A). The oscillations were re- 
ported as early as the 19th century and are evident for 
instance in the extremely high catches by fisheries every 
fourth year-- . This both economically and conservation- 
ally important phenomenon has been ascribed either to 
transient effects or to stochastic infiuences^^, to depen- 
satory predatiorfl^, to fishin^i^, or to genetic effect J^^l^^, 
but all of these explanations are still not fully convincing 



The sockeye salmon return to spawn in their native 
stream or lake in late summer and then die, which means 
that the salmon generations do not overlap. The hatched 
fry migrate downstream in the following spring and feed 
for one season in large freshwater lakes, before they mi- 
grate to the ocean, where they spend the next two and a 
half years. The carcasses of the adult salmon are decom- 
posed and provide an important phosphorus input into the 
rearing lakes of the fry. Since the oscillations of sockeye 
salmon originating from different lakes are not in phase, 
we have clear evidence that the relevant processes caus- 
ing the phenomenon of cyclic dominance take place in the 
rearing lakes rather than in the ocean. 
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Figure 1: Time series of sockeye salmon abundance. 

Total escapement (number of male and female adults that es- 
cape the fisheries) of the late Shuswap Lake run. Shuswap lake 
is one of the most important rearing lakes of juvenile sockeye 
salmon. Data for 5 other lakes are provided in Appendix A. 
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2 Model and results of computer 
simulations 

It is our aim to develop a generic model capturing only the 
essential mechanisms required for the occurrence of cyclic 
dominance. As such, the model is kept very simple and is 
not designed to quantitatively predict the population dy- 
namics of all species in the corresponding lake ecosystems. 

The model uses standard continuous population dynam- 
ics equations for sockeye fry, Sn(i), their predator (e.g. 
rainbow trout), Pn{t), and their zooplankton food (mainly 
daphnia), z„(t), during the growth season from spring 
{t = 0) to fall {t = T) in year n, combined with a rule 
for calculating the three population sizes at the beginning 
of the next season as a function of the population sizes 
at the end of the previous season(s) (see Materials and 
Methods). The sockeye fry population at the end of a sea- 
son gives rise to the number of spawning adults 3 years 
later, which in turn determines the number of sockeye fry 
in the following spring, and represents a nutrient input 
that affects the carrying capacity of the zooplankton. A 
small fraction of the salmon stays in the ocean for one 
more year and mature at age 5, thus causing a mixing 
between the four brood lines of a spawning population. 
An even smaller fraction matures at age 3, but since these 
fish are predominantly small-sized male^ that do not in- 
fluence the number of fertilised eggs, they are neglected in 
this study. 

According to empirical observations ~, the zooplank- 
ton level at the end of one year has no effect on the 
following year. However, its carrying capacity is a func- 
tion of the nutrients provided by the adult salmon (and 
thus a function of the initial fry biomass of that year): 
Kn = Kq + /(s„(0)), where Kq is the carrying capacity in 
the absence of sockeye spawners, and / is a function that 
increases with the number of spawners. The zooplank- 
ton is initialised each year with its carrying capacity. The 
predator biomass does not change over the winter season 
in our model. 

Fig. 2a shows a time series of the biomass of the sockeye 
fry at the end of the growth season. The first 300 years 
are cut off to show only the stable periodic oscillation with 
one strong year followed by one intermediate year and two 
weak years, just as in the empirical data of most sockeye 
stocks showing cyclic dominance (Fig. 1 and Appendix A) . 
When the simulation parameters are chosen differently, 
the difference between the strong and the intermediate 
year may become less pronounced, or the order may even 
become reversed. This is also observed in the empirical 
data of some stocks. 

When a parameter is varied, for instance the constant 
fraction of the zooplankton carrying capacity, Kq, the dy- 
namical pattern may change. Fig. 3 shows the biomass 
s„ (T) of the sockeye fry at the end of the season for dif- 
ferent values of Kq, from year t =1000 to year 1100. For 
small Kq, all sockeye lines are equally strong. This means 
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Figure 2: Simulated time series. Stable attractor of sock- 
eye salmon dynamics generated by computation of the time- 
continuous model, showing oscillations with period 4 and the 
same pattern in the sequence of salmon abundance as in the 
empirical data. Only the biomasses at the end of each growth 
season are shown, a: deterministic model, b: dynamics with 
up to 50 percent fluctuations in the survivability of the sock- 
eye salmon in the ocean (see Discussion) . Parameter values for 
these simulations are as explained in Materials and Methods. 



that the dynamics reaches a fixed point and that there 
is no cyclic dominance. With increasing Kq, the fixed 
point eventually becomes unstable. First a bifurcation to 
quasiperiodic behaviour occurs, and then the frequency of 
the oscillation becomes locked at 4. 

3 Linearised Theory 

In order to understand and interpret these observations, 
we first note that the continuous population dynamics dur- 
ing the season, together with the matching conditions ap- 
plied between two seasons, can be viewed as a discrete 
map (a so-called Poincare map), giving the sockeye and 
predator biomasses at the end of one year as function of 
the biomasses at the end of the previous years. 

To obtain this map, one first has to integrate the dy- 
namical equations over one season, giving s„(T) and Pn(T) 
as a function of s„(0) and Pn(0). The zooplankton can be 
eliminated since the initial condition of the zooplankton 
depends on Sn(0) only. Next, one expresses Sn-\-\ (0) and 
Pn+i{^) as a function of s„_3(T), s„_4(T) and pn(T), us- 
ing the matching conditions. 

The mechanism which generates the population oscilla- 
tion is based on the nature of the instability of the sta- 
tionary state s„(r) = s*{T), pn{T) = p*{T) of our sys- 
tem. The corresponding bifurcation can be investigated 
in terms of a linear stability analysis. Close to the fixed 
point, the dynamics can be approximated by linear terms. 
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Figure 3: Bifurcation diagram of the biomass of the 
sockeye fry at the end of the season. Bifurcation parame- 
ter is the constant fraction of the zooplankton carrying capac- 
ity Kq. For each value of Ko, 100 consecutive data points are 
plotted. Whenever only a small number of different points is 
visible for a value of the bifurcation parameter, this indicates 
periodic oscillations or a fixed point. 



Denoting the distance of the biomasses from their fixed 
point value by Ssn — s„(T) — s*(T) etc, the linear approx- 
imation of this map has the form 
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with positive parameters ms,mp,a,b,c,ei. The general 
structure of the matrix in Eq. 1 is determined only by the 
matching conditions and is valid for any time-continuous 
model applied during the seasons. The equations of mo- 
tion of the latter only determine the numerical values of 
the parameters in Eq. (1). 

The first line of the matrix describes the influence on 
the sockeye fry of year n -I- 1 of the sockeye fry of year 
(n — 3) and (n — 4) (which are the parents of the fry in 
year (n -I- 1)), and of the predator population, and ei 
are positive, since more parents imply more offspring, —a 
is negative, since more predators imply less fry. The other 
nontrivial line of this matrix, the last line, describes the 
influence on the predators of year (n + 1) of the sockeye 
fry of year (n — 3) and (n — 4) (which are the parents of 
the predator's food), and of the predator population in 
the previous year. All three parameters are positive, since 
more food implies more predator growth and since more 
predators in one year give rise to more predators in the 
next year. 

The eigenvalues of this matrix determine the nature of 
the dynamics near the bifurcation. When all eigenvalues 
have an absolute value smaller than 1, the fixed point is 
stable, and the dynamics converges to this fixed point. 
When the absolute value of one or more eigenvalues is 
larger than 1, the fixed point is unstable, and the dynamics 



approaches a different attractor. 

In order to understand the dynamics resulting from this 
matrix, we first consider the case that the matrix element 
ei and the product a{b + c) vanish. (Note that since all 
parameters are positive, the latter requires both products 
ab and ac to vanish.) This means that all salmon return at 
age 4, and that the trout have a good choice of other food 
so that their dynamics is independent of that of the salmon 
fry. In this case the eigenvalues of the matrix are mp, 0, 
and the four fourth roots of iris. Since the four salmon 
lines are independent from each other in this case, the 
sequence Sst has trivially the period 4 and simply iterates 
the initial four values, with an amplitude decreasing for 
rus < 1 and increasing otherwise (and with the trout being 
independent of the salmon). When rris is increased from 
a value smaller than 1 to a value larger than 1, all four 
eigenvalues ml^^ cross the unit circle simultaneously, and 
the fixed point becomes unstable. This degeneracy is lifted 
when the parameters ei and a{b + c) are made nonzero. 
As long as these parameters are not large, one can expect 
the four main eigenvalues to remain close to the real and 
imaginary axis, respectively, implying a (possibly damped) 
oscillation with a period close to 4. 

The type of bifurcation that occurs when the fixed point 
becomes unstable depends on which eigenvalue crosses 
first the unit circle as a control parameter is increased. 
In order to determine the type of the bifurcation, we eval- 
uated the parameters of the matrix numerically from our 
computer simulations of the time continuous model near 
the bifurcation, and we calculated the eigenvalues A of the 
matrix. Figure 2] shows that the bifurcation is caused by 
a pair of complex conjugate eigenvalues crossing the unit 
circle, indicating a Neimark Sacker bifurcation (the dis- 
crete version of the Hopf bifurcation) . Since the dominant 
eigenvalues are close to ±i, the period of the resulting os- 
cillation is close to 4. The bifurcation is mainly driven 
by the parameters tUs and, to a minor extent, a, which 
correspond to reproduction of the 4-year old sockeyes and 
predation by the rainbow fronts, respectively. The fixed 
point becomes unstable when either of these parameters 
increases. 

Now, it is known from the theory of bifurcations that 
if the period at such a bifurcation is close to 4, there oc- 
curs a strong resonance, which means that the period be- 
comes locked exactly at the value 4 not far beyond the 
bifurcation point. In contrast to conventional frequency 
locking, a strong resonance is due to nonlinearities that 
are of the same order as the leading nonlinearity, and fre- 
quency locking therefore occurs over a much wider range 
of parameter^i^. 

The two complex conjugate eigenvalues are the first ones 
to cross the unit circle when a and ei are small non- 
negative numbers and when nip is well below 1. Since ei 
and a{b + c) are small compared to m^, the period of the 
resulting oscillation remains close to 4, and the locking at 
period 4 therefore occurs not far beyond the bifurcation. 
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Figure 4: Eigenvalues of the matrix from Eq. 4, as 
obtained from a linear stability analysis of the time- 
continuous model. Ko increases from 6.4 to 7.5 as the eigen- 
values move outwards. The zoom shows that the eigenvalues 
close to the imaginary ajds are the first ones to cross the unit 
cycle. 



When these conditions are not met, the unit circle is typ- 
ically first crossed at -1, and a period-doubling (or flip) 
bifurcation occurs. 



4 Discussion 

All these considerations lead to two basic conditions un- 
der which strong resonance can be observed: First, an 
increase in the number of spawning salmon must lead to a 
sufficiently strong increase in this number four years later 
(i.e. rUs must be large enough). If the number of fry mi- 
grating to the ocean is dominated by other factors, such as 
a strongly limited carrying capacity for the fry, the popu- 
lation will be at a fixed point rather than on the oscillating 
side of the bifurcation. Second, the four salmon lines must 
be coupled in order for the Neimark Sacker bifurcation to 
occur, rather than a period doubling bifurcation. In our 
model, this coupling is due to a fraction of sockeye re- 
turning at age 5 instead of age 4, and, more importantly, 
due to the predator being sufficiently strongly coupled to 
sockeye dynamics. 

Previous studies of salmon dynamics, used in salmon 
management, are based on the Ricker modeP"^ or the 
Larkin modeP^'^^'^. Because both models have a 
strongly limited carrying capacity for the fry, and because 
they do not include explicitly the predator dynamics, the 
only bifurcation occurring in those models is the flip bi- 
furcation. Building on both modelling approaches, Myers 
et have demonstrated that an unstable (decaying) 

period-4 oscillation can be excited by stochastic driving. 
However, it remains unclear over how long the oscilla- 
tions can be maintained in that study. Furthermore, the 
frequently observed sequence of population numbers de- 
scribed in section 2 (and stably seen in our simulations. 
Fig. only occurs episodically in the time series shown 



due to the stochastic nature of the approach. 

In principle, the coupling between the four salmon lines 
can also occur through the food of the sockeye fry. How- 
ever, there is no empirical evidence that there is a negative 
effect of a strong sockeye year on the daphnia population 
in the following year that limits the growth of the sockeye 
fry in that yea.^^. 

The two conditions for a strong resonance of period 4 
fit together with the empirical observation that it occurs 
only in large oligotrophic lakes, such as those of the Fraser 
River basin. The smaller and ultra-oligotrophic lakes in 
the outer coast regions of British Columbia do not show 
these oscillations. Their nutrient content is neither large 
enough to raise large fry numbers, nor to allow for salmon 
predators to become strong. 

On the other hand, cyclic dominance cannot be ex- 
pected in nutrient-rich lakes, because the spawning adults 
would not be an important nutrient input, and because 
the coupling to the predator would not be strong enough 
in a situation where there were more predator species and 
more prey species for these predators. 

Since there are large fluctuations in the proportion of fry 
that survive to return to their nursery lake, we included 
noise in the matching condition for the sockeyes (first Eq. 
3 in Materials and Methods) in order to determine with 
how much noise superimposed on the deterministic dy- 
namics the period-4 oscillation can persist. Fig. 2b shows 
a data series generated with 50 percent noise in ocean 
survivability superimposed to the deterministic dynamics. 
The oscillation is still clearly visible, although the sys- 
tem has a phase shift every 300 years on average. With 
less noise, the phase shift occurs less often, and with 100 
percent noise, the cyclic dominance vanishes in our simu- 
lations. When a large perturbation acts only for a limited 
time, the oscillation quickly recovers afterwards. In fact, 
the recovery following the blocking of the Fraser River 
migration routes early in the last century, can be seen 
in the non-dominant brood lines of the Shuswap stock in 
Fig. 1 and in Fig. Al in Appendix A. Unfortunately, the 
year 2009 has seen another large perturbation, with most 
Fraser sockeye expected to return that year not surviving 
in the ocean, so that the expected strong escapement did 
not occur. 

The age composition of sockeye stocks was also found 
to strongly infiuence the potential of the system to show 
a strong resonance. The fraction of adult sockeyes that 
returns at the age of 5 years instead of 4 years, e, was 
set to 0.1 in the simulations and we could show that the 
resonance appears for < e < 0.3, but for very small val- 
ues of e (< 0.02) the non-dominant lines disappear com- 
pletely with increasing Kq. The resonance condition was 
best met at e « 0.2, where the resonance occurred nearly 
directly after the Neimark Sacker bifurcation. The param- 
eter e is also measured in the real population^^Ql^ Yot the 
sockeye salmon populations of the Fraser River it is ap- 
proximately 0.08, while for the less productive outer coast 
lakes of British Columbia (where cyclic dominance is not 
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observed) it is between 0.56 and 0.76. 

Some of the sockeye populations of the Bristol Bay area 
(Alaska), most notably the Kvichak River stocks, also 
exhibit strong oscillation^^illlll despite a broader distri- 
bution of the age at spawning than in the Fraser River 
stock^2llll!. However, the oscillations are not as regular, 
with maxima of the populations occurring every fourth or 
fifth year. This is consistent with our model, as it indicates 
a quasiperiodic oscillation rather than a fixed oscillation 
period associated with a strong resonance. 

Our results do not rule out additional mechanisms such 
as depensatory fishing^'^ or genetic effects ^i^, which would 
reduce the population sizes of weak lines to values smaller 
than those resulting from our model. However, these ad- 
ditional assumptions are not needed to explain the occur- 
rence of cyclic dominance in the first place. 



5 Materials and Methods 
5.1 Model equations 

The dynamics of the biomass of sockeye fry s„(t), of their 
predator p„(i), and of their zooplankton food Zn{t) in year 
number n during the growth season from spring {t = 0) to 
fall {t = T) are given by the following equations: 



on s„+i(0)). 



—a 



Snjt) -Pnit) 



ds ■ Sn{t) 



Zn{t) = Zn{t) • 1 - 



dt 



Zn{t) ■ Sn{t) 



(2) 



l+Cs ■ Sn{t) + Zn{t) 



Sn{t) ■Pn{t) 



Pnit) = \-aps— T— —--dp -Pnit) 

1 + Cp - Pnit) + Sn(t) 



dt 



The feeding terms include saturation at high prey densi- 
ties, and a predator interference term in the denominator 
(Beddington functional respons^^H HH) _ Similar predator 
interference terms are also used by other modeller^^H 

The matching conditions used to determine the 
biomasses of the species at the beginning of the next sea- 
son from their values at the end of the previous season(s) 
are given by 



s„+i(0) = 7((l-e)s„_3(T) + es„_4(T)) 

Z„+i(0) Kn+l 
Pn+l(0) = Pn{T), 



(3) 



with e the proportion of surviving sockeye that return to 
their native lakes at the age of 5 to spawn and die. The 
carrying capacity of the zooplankton in the next season, 
Kn+i, shows a saturating dependence on the nutrient in- 
put due to the number of spawning adults (and thereby 



s«+i(0) 
5o + s„+i(0) 



(4) 



5.2 Parameters used for the computer 
simulations 

A — 0.85 denotes the assimilation efhciency of ingested 
prey biomass of carnivores. Ugf — 10 and aps = 1.6 are 
the maximal per unit biomass ingestion rates of salmon 
and predators, respectively, dg — 1 and dp = 0.16 repre- 
sent biomass loss due to respiration and mortality. The 
metabolic rates of the predators are smaller than the re- 
spective rates of the salmon to account for the larger body 
size of the predators ^^Sl; we assume the predators to be ap- 
proximately 1500 times heavier than the sockeye fry. The 
interference parameters are set to Cg = 1.0 and Cp — 0.2. 
The parameter 7 in Eq. 3 summarises survival from smolt 
to adult fish (including ocean survivability and loss to fish- 
eries), spawning success, and egg to fry survival. The vari- 
ous factors are estimated from empirical data'^^'^^'^lland 
yield 7 = 0.4. The 3 parameters determining the carry- 
ing capacity Kn of the zooplankton in year n are set to 
Ko = 15, (5 = 10, and Sq = 5. 



Acknowledgements 

This work was supported by the German Research 
Foundation (DFG) under contract number Dr300/7 and 
Br2315/9-l. Kim Hyatt, Jeremy Hume, and Carl Wal- 
ters commented on an earlier version of this manuscript. 
Empirical data were kindly provided by Tracy Cone (es- 
capement data) and by Kim Hyatt (age composition of the 
sockeye stocks). This collaboration started at the Com- 
plex Systems Workshop in Fairbanks, Alaska, 2007. The 
2008 Cultus Lake Modelling Workshop, which had an im- 
portant impact on this project, was hosted by Fisheries 
and Oceans Canada. 



References 

[1] Ehon C, Nicholson M (1942) The ten-year cycle in 
numbers of lynx in Canada. J. Anim. Ecol. 11:215- 
244. 

[2] Odum EP (1953) Fundamentals of Ecology (W. B. 
Saunders, Philadelphia). 

[3] Blasius B, Huppert A, Stone L (1999) Complex dy- 
namics and phase synchronization in spatially ex- 
tended ecological systems. Nature 399:354-359. 

[4] Hanski I, Turchin P, Korpimaki E, Henttonen H 
(1993) Population oscillations of boreal rodents - 
Regulation by mustelid predators leads to chaos. Na- 
ture 364:232-235. 



5 



[5] Ludwig D, Jones DD, HoUing CS (1978) Qualitative 
analysis of insect outbreak systems - spruce budworm 
and forest. J. Anim. Ecol. 47:315-332. 

[6] Royama T (1984) Population dynamics of the 
spruce budworm choristoneura-fumiferana. Ecol. 
Mon. 54:429-462. 

[7] Ricker WE (1950) Cycle dominance among the Eraser 
sockcye. Ecology 31:6-26. 

[8] Townsend CR (1989) Population cycles in freshwater 
fish. J. Fish. Biol. 35 (Supplement A):125-131. 

[9] Ricker WE (1997) Cycles of abundance among Eraser 
River sockeye salmon (Oncorhynchus nerka). Can. J. 
Fish. Aquat. Sci. 54:950-968. 

[10] Rounsefell GA, Kelez GB (1938) The salmon and 
salmon fisheries of Swiftsure Bank, Puget Sound, and 
the Eraser River. Bull. Bur. Fish. 49:692-823. 

[11] Myers RA, Mertz G, Bridson JM, Bradford MJ 

(1998) Simple dynamics underlie sockeye salmon (On- 
corhynchus nerka) cycles. Can. J. Fish. Aquat. Sci. 
55:2355-2364. 

[12] Larkin PA (1971) Simulation studies of the Adams 
River sockeye salmon, Oncorhynchus nerka. J. Fish. 
Res. Board Can. 18:1493-1502. 

[13] Walters CJ, Staley MJ (1987) Evidence against the 
existence of cyclic dominance in Eraser River sockeye 
salmon (Oncorhynchus nerka). In: Sockeye salmon 
(Oncorhynchus nerka) population biology and fu- 
ture management. Can. Spec. Publ. Fish. Aquat. Sci. 
96:375-384. 

[14] Levy DA, Wood CC (1992) Review of proposed mech- 
anisms for sockeye salmon population cycles in the 
Eraser River. Bull. Math. Biol. 54:241-261. 

[15] Walters CJ, Woodey, JC (1992) Genetic models for 
cyclic dominance in sockeye salmon (Oncorhynchus 
nerka). Can. J. Fish. Aquat. Sci. 49:281-292. 

[16] Kuznetsov YA (2004) Elements of Applied Bifurca- 
tion Theory (Springer- Verlag, New York). 

[17] Hume JMB, Shortreed KS, Morton KE (1996) Ju- 
venile sockeye rearing capacities of three lakes in 
the Eraser River system. Can. J. Fish. Aquat. Sci. 
53:719-733. 

[18] Martell SJD, Walters CJ, Hilborn R (2009) Retro- 
spective analysis of harvest management performance 
for Bristol Bay and Eraser River sockeye salmon. Can. 
J. Fish. Aquat. Sci. 65:409-424. 

[19] Marsden AD, Martell SJD, Sumaila UR (2009) Retro- 
spective bioeconomic analysis of Eraser River sockeye 
salmon fishery management. Fish. Res. 95:32-41. 



[20] Healey M (1987) The adaptive significance of age 
and size at maturity in female sockeye salmon (On- 
corhynchus nerka) . p. 110-117 in: Smith HD, Margolis 
L, Wood CC (Eds.), Sockeye Salmon (Oncorhynchus 
nerka) Population Biology and Euture Management. 
Can. Spec. Publ. Fish. Aquat. Sci. 96, Eisheries and 
Oceans Canada, Ottawa. 

[21] Eair LE (2003) Critical elements of Kvichak River 
sockeye salmon management. Alaska Eish. Res. Bull. 
10: 95103. 

[22] Rogers LA, Schindler DE (2008) Asynchrony in pop- 
ulation dynamics of sockeye salmon in southwest 
Alaska. Oikos 117:1578-1586. 

[23] West EW, Eair LE (2006) Abundance, age, sex, and 
size statistics for pacific salmon in Bristol Bay, 2003. 
Fishery Data Series No. 06-47, Alaska Department of 
Eish and Game, Division of Sport Eish and Commer- 
cial Eisheries, Anchorage. 

[24] West EW, Eair LE, Baker T, Morstad S, Weiland K, 
Sands T, Westings C (2009) Abundance, age, sex, 
and size statistics for pacific salmon in Bristol Bay, 
2004. Fishery Data Series No. 09-51, Alaska Depart- 
ment of Eish and Game, Division of Sport Fish and 
Commercial Eisheries, Anchorage. 

[25] Beddington JR (1975) Mutual interference between 
parasites or predators and its effects on searching ef- 
ficiency. J. Anim. Ecol. 44:331-340. 

[26] Skalski GT, Gilham JE (2001) Functional responses 
with predator interference: viable alternatives to the 
HoUing type H model. Ecology 82:3083-3092. 

[27] Walters CJ, Christensen V (2007) Adding realism to 
foraging arena predictions of trophic flows in Ecosim 
ecosystem models: shared foraging arenas and bout 
feeding. Ecol. Mod. 209:342-350. 

[28] Brown JH, Gillooly JE, Allen AP, Savage VM, West 
GB (2004) Toward a metabolic theory of ecology. 
Ecology 85:1771 1789. 

[29] Pauley GB, Risher R, Thomas GL (1989) Species pro- 
files: life histories and environmental requirements of 
coastal fishes and invertebrates (Pacific Northwest) 
- Sockeye salmon. US Dep. Int., Fish and Wildlife 
Service, Biological Report 82(11.116). 



6 



A Supporting empirical data 



Escapement data for sockeye stocks showing cyclic dominance Many spawning stocks of sockeye salmon in 

the upper Frascr River basin are not stationary but have been dramatically increasing in size over the last decades (e.g. 
Fig. SI a and b), following a massive disturbance of the stocks at the beginning of the 20th century. Nevertheless, 
the oscillatory pattern has clearly emerged again. 




1940 1950 I960 1970 1980 1990 2000 
Year 



Figure Al: Escapement data of for sockeye spawing stocks in the Eraser River basin, a: Quesnel Lake, b: Stuart Lake, c: 
Bowron Lake, d: Seton Lake, e: Chilko Lake. 
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